Robust dynamical decoupling for quantum computing and quantum memory 
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Dynamical decoupling (DD) is a popular technique for protecting qubits from the environment. 
However, unless special care is taken, experimental errors in the control pulses used in this technique 
can destroy the quantum information instead of preserving it. Here, we investigate techniques for 
making DD sequences robust against different types of experimental errors while retaining good 
decoupling efficiency in a fluctuating environment. We present experimental data from solid-state 
nuclear spin qubits and introduce a new DD sequence that is suitable for quantum computing and 
quantum memory. 
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An obstacle against hiffh-precision quantum control is 
the decoherence process [l[ . The ability to preserve quan- 
tum behavior in the presence of noise is essential for the 
performance of quantum devices, such as quantum mem- 
ories 0, where one wishes to store a quantum state, and 
quantum computers, where the quantum information is 
processed Q. 

A promising strategy developed to avoid decoherence 
is the dynamical decoupling (DD) method which aims 
to reduce decoherence times by attenuating the system- 
environment interaction. Since DD does not require aux- 
iliary qubits or measurements, it can be used as an eco- 
nomical alternative to complement quantum error cor- 
recting codes . Decoupling schemes were originally de- 
veloped in the framework of nuclear magnetic resonance 
(NMR) 0. In DD, a sequence of control fields is pe- 
riodically applied to a system in cycles of period Tc, in 
order to refocus the system-environment evolution. The 
delay between pulses is one of the relevant parameters 
of a sequence When Td is shorter than the correla- 
tion time Te of the environment, the preservation of a 
single qubit state is ideally possible even for the most 
general dephasing environment The effectiveness of 
a decoupling scheme depends crucially on the repetition 
rate with which the pulses can be applied. Recent ex- 
periments have successfully implemented DD methods 
and demonstrated the resulting increase of the coherence 
times 0,0 However, systems with fast fluctuating 
environment, with short Te, put enormous demands on 
the hardware. They may be encountered in a wide range 
of quantum information processing implementations [11 1 
and represent the most challenging regime for avoiding 
decoherence. 

Apart from preserving the state of a quantum memory 
as long as possible, DD can also be used to keep a quan- 
tum state coherent while logical operations are performed 
on it. In this case, one needs to refocus the system- 
environment interactions while not refocusing the desired 
system-system interactions needed for multi-qubit gates, 
as discussed in Ref. [l^l- Suitable DD sequences for 
quantum computing should therefore keep the number of 



refocusing pulses small to allow the computational opera- 
tions and limit power deposition on the sample. Since the 
precision of any real operation is finite, the performance 
of experimentally accessible DD sequences is limited also 
by the pulse errors 0, E, El- This effect must be min- 
imized by choosing a decoupling scheme that is robust 
against pulse imperfections. 

As the spacing between the pulses is reduced, the refo- 
cusing gets more effective, since the state of the environ- 
ment appears more static on the time scale of the pulse 
spacing. Ideally, this can be extended indefinitely, until 
the system decouples completely from the environment 
in the limit of infinitely short pulse spacing. In practice, 
however, the observed decay rate goes through a mini- 
mum and subsequently increases again. In this regime, 
the signal is destroyed primarily by imperfections of the 
refocusing pulses. 

We first consider this limiting case, where the effect 
of the environment is negligible and we need to discuss 
only the effect of control errors. The dominant cause of 
these errors is, in most cases, a deviation between the 
actual and the ideal amplitude of the control field. The 
result of this amplitude error is that the rotation angle 
deviates from tf, typically by a few percent. A possible 
approach for compensating these errors is the use of com- 
posite pulses [l3l , which generate rotations that are close 
to the target value even in the presence of amplitude er- 
rors. In this case, the error correction is done "inside" the 
pulse. Alternatively, it is possible to design the sequence 
in such a way that the error introduced by one pulse is 
compensated by subsequent pulses. We refer to the for- 
mer approach as using robust pulses and to the second 
approach as using robust or self-correcting sequences. 

Apart from reducing the coherence time in a certain 
parameter range, pulse imperfections also make some DD 
sequences asymmetric with respect to the initial con- 
ditions. We illustrate this in Fig. [TJ where we plot 
the signal of the carbon nuclear spin magnetization in 
the Adamantanc molecule measured for two different 
DD sequences and initial conditions. Considering first 
the Carr-Purcell-Mciboom-Gill (CPMG) sequence, which 
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Fig. [31 causes a slower decay for the same parameters, 
indicating a partial error compensation over the cycle. 
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Figure 1: Normalized magnetization obtained experimentally 
for two basic NMR sequences: one non-robust against errors, 
CPMG, and one self-correcting sequence, XY-4. 



consists of identical 7r-pulses. It was originally desi gned 
to preserve a single component of the magnetization [1^ . 
Here, the decoherence time is sa 100 ms if the Bloch vec- 
tor of the qubit is initially oriented parallel to the pulse 
axis (longitudinal state). Under these conditions, pulse 
errors do not affect the coherence llSl - they even pro- 
vide additional stabilization 0, [l6l. Il7|. In contrast, if 
the initial condition is perpendicular to the direction of 
the pulses (transverse state; blue squares in Fig. [Ij, the 
errors of the individual pulses accumulate and lead to 
a rapid decay: the signal completely vanishes after « 50 
pulses (« 5 ms). A similar behavior is found for the UDD 
sequence which also uses rotations around a single 
axis [ill. 

The second DD sequence represented in Fig. [T] is the 
XY-4 sequence which consists of tt pulses applied 
along the x and y axes It performs much more symmet- 
rically with respect to the initial condition: both initial 
states decay on a time-scale of « 10 ms. The observed 
time scale also shows that this sequence has a built-in 
partial error compensation. 

For preserving an unknown quantum state, the appro- 
priate performance measure should not depend on the 
initial condition. A common choice for quantifying the 
performance of quantum operations is then the fidelity 



F 



\Tr{AB'<)\ 



^Tr{AA'i)Tr{BB^) 



(1) 



Here, A is the target propagator (unity for the examples 
discussed here) and B is the actual propagator generated 
by the real pulse sequence. In Fig. [2l we numerically 
simulate the fidelity decay during the two DD sequences 
discussed above due only to errors in the fiip angles of the 
pulses, neglecting environmental effects. For the CPMG 
sequence (blue circles), we observe a rapid decay to the 
limiting fidelity of « 0.65. At this point, the system 
reaches the completely disordered state p oc 1, with the 
order dephased in the inhomogeneous field distribution. 
The fast decay for CPMG is experimentally manifested 
by the large asymmetry shown in Figure [T] The XY- 
4 sequence, however, represented by the red squares in 
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Figure 2: Simulation of fidelity decay due to pulse errors. 
The blue circles represent the decay of the fidelity due only to 
pulse errors for the CPMG sequence and the red squares for 
the XY-4 sequence. The fidelity is calculated as an average 
over a Gaussian distribution of flip angles; the width was 10% 
of the nominal flip angle. 

In order to correct errors of non-robust sequences or 
even to improve the error tolerance of self-correcting se- 
quences, it is possible to replace the simple tt pulses by 
composite pulses [31 that are more robust against errors. 
We tested different class- A composite 7r-pulses, which 
produce compensated rotations for any initial condition, 
and found that the pulse used in Rcf. [l3| 

(7I")7r/6+</. - (tt),/, - (7r)^/2+0 - (tt)^ - (7r)^/6+07 (2) 

to which we refer here as the Knill pulse, is the most 
robust against flip-angle errors and off-resonance errors, 
which are the leading errors in many experimental situa- 
tions. The Knill pulse is equivalent to a robust tt rotation 
around the axis defined by (p followed by a —■k/'Ii rotation 
around the z axis. For cyclic sequences, which always 
consist of even numbers of tt rotations, the effect of the 
additional z rotation vanishes if the flip angle errors are 
sufficiently low. We inserted these composite pulses into 
different DD sequences to test their performance. 

Since the Knill pulse consists only of tt rotations, it can 
also be used as a DD sequence: instead of concatenating 
the pulses directly, we inserted delays between them and 
obtained thus a new DD sequence with the rotation axes 
given by Eq. 0: KDD^ = t/2 - {tt)^/q+^ - t - (tf)^ - 
T - (7i')7r/2+0 - T - (7r)0 - T - (7r)^/g+^ - t/2. To furthcr 
improve the robustness of the sequence, we also extended 
it by combining 5-pulsc blocks shifted in phase by 7r/2, 



such as [KDD0 - KDD^+7r/2] , where the lower index 
gives the overall phase of the block. We will refer to the 
cyclic repetition of these 20 pulses as the KDD sequence. 

Fig. [3] summarizes how sensitive the different DD se- 
quences are to the two main pulse imperfections (no de- 
phasing environment): Each panel shows the simulated 
fidelity F after 100 pulses with a combination of fre- 
quency (vertical axis) and amplitude errors (horizontal 
axis). Fidelities > 0.95 arc color coded. For CPMG 
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Figure 3: Error tolerance of difTerent DD sequences. The up- 
per row shows the calculated fidelity F for standard DD se- 
quences, while the lower row shows the results for the CPMG 
and PDD=XY-4 sequences when the tt pulses are replaced 
by Knill pulses. The last panel corresponds to the KDD se- 
quence, also based on the Knill pulse. Each panel shows the 
fidelity after 100 pulses as a function of flip angle error and 
offset errors. The regions where the fidelity is lower than 0.95 
are shown in white. The highest contour level is _F = 0.999. 



(first panel), the fidelity after 100 pulses drops to <0.95 
even for very small flip angle errors or offsets. The 
concatenated sequences (CDD-n) were constructed as 
[C'DDn^i -X- CDDn-i - YY and yCDDn-i ~X~ 
CDDn~i —Y— sjCDDn^iY' for the standard and sym- 
metrized versions respectively. For n = 1, we have CDDi 
= XY-4 = FDD. The basic sequence XY-4 is defined as 
[Td — X — Td — YY for the standard case. The symmetrized 
sequence is [Td/2 - X - Td - Y - Td/2]^. The 2"'* and 
3'"'^ panel in Fig. [3] show the corresponding results for 
the standard (asymmetric) XY-4 and CDD2 sequence. 
Clearly, they are much less susceptible to flip angle er- 
rors than the CPMG sequence. As shown in the lower 
panels, a further significant improvement is achieved if 
the TT-rotations are replaced by Knill pulses, at the ex- 
pense of increased power deposition (by a factor of 5). 
The best performance, with similar power deposition as 
for the panels in the upper row, is achieved with the KDD 
sequence. 

Now we start looking at the effect of imperfect DD 
pulses on a fluctuating environment. We tested the se- 
quences discussed above in an experimental setting. As 
the system qubit, we used ^'^C nuclear spins in the CH2 
groups of a polycrystalline Adamantane sample. The 
natural abundance carbon spins are surrounded by 
nuclear spins acting as a rapidly fluctuating environment 
The bath correlation time was ~ 100/j,s, the pulse 
length was 10.6/xs and the delays between the pulses were 
varied from 5/Lts to 150^s. Under our conditions, the in- 
teraction between the carbon nuclei can be neglected and 
the decoherence mechanism is a pure dephasing process 
The experiments were performed on a home-built 
300 MHz solid-state NMR spectrometer. 

Figures 2] and [5] summarize the experimental results 



by plotting the relaxation times as a function of the 
duty cycle (total irradiation time divided by total time) . 
The relaxation times are defined as the 1/e decay time 
of the magnetization. In the first set of experiments. 
Fig. UJ we compare different members of the XY fam- 
ily. The sequences were constructed, as explained in [2^ , 
using the standard definition of XY-4 and the its sym- 
metric definition, introduced above. For all sequences, 
the symmetric version performs significantly better, with 
relaxation times ~ 5 times longer than for the asym- 
metric version, irrespective of the initial condition. This 
can be attributed to the fact that in the symmetric se- 
quences, all odd-order terms of the Magnus expansion 



vanish [21|. This is a significant advantage, considering 
that the power deposition and the complexity of the se- 
quences are identical. 

In the second set of experiments. Fig. [5l we compare 
standard CDD sequences against CDD sequences with ro- 
bust pulses and symmetric timing. For low duty cycles, 
standard CDD sequences perform better. This is due to 
the shorter cycle time of the standard sequence if con- 
stant duty cycles are compared. Using non-robust pulses 
may therefore be the preferred option if DD sequences 
are applied parallel to gate operations. However their 
performance always saturates or decreases with increas- 
ing duty cycle under the present experimental conditions, 
while the performance of sequences with robust pulses 
continuous to improve. Thus, if only the preservation of 
a quantum state is required, without considering power 
deposition, the sequences with robust pulses provide the 
best performance. The performance of the sequences 
with robust pulses shows no significant dependence on 
the concatenation level. Apparently, the additional cor- 
rections due to the concatenation are not required in this 
case. 
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Figure 4: Comparison of the decoupling performance of sym- 
metric vs. asymmetric XY sequences. 



For a large range of duty cycles, we find that KDD pro- 
vides the best performance. At low duty cycles the per- 
formance of KDD is comparable to that of self-correcting 
sequences without robust pulses, indicating that the er- 
rors of the individual pulses compensate over a cycle. For 
high duty cycles, instead of saturating, the relaxation 
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Figure 5: Comparison of the decoupling performance of stan- 
dard DD sequences (FDD, CDD2) with the same sequences 
using composite pulses and the self-correcting sequences XY- 
16 and KDD. 

time continues to increase, as in the case of sequences 
with robust pulses. This shows that KDD has suitable 
properties for computing and memory applications. 

In summary, we have considered the problem of pro- 
tecting a qubit in the presence a fast fluctuating spin- 
bath and pulse imperfections. Different strategies for 
correcting pulse errors were investigated and verified ex- 
perimentally in DD experiments. We observed that the 
symmetrization of sequences is an important feature since 
this leaves the power deposition and the complexity of 
the sequence constant but always leads to better perfor- 
mance. The best sequences that arc suitable for parallel 
application of quantum gate operations are the symmet- 
ric self-correcting sequences. However, their performance 
saturates at higher duty cycle, while the performance of 
sequences with robust pulses continues to improve un- 
der our experimental conditions. Thus, if the objective is 
only to preserve a quantum state, the best performance 
is obtained at high duty cycles, using robust pulses. We 
also introduced a new DD sequence, which combines the 
useful properties of robust sequences with those of robust 
pulses and can thus be used for both quantum comput- 
ing and state preservation. This new sequence contains 
rotations around different axes in the xy plane, not only 
around the more conventional directions x and y. We be- 
lieve that mixing nontrivial directions will be a helpful in- 
gredient for developing new robust dynamical decoupling 
sequences for future quantum information applications. 
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